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1.  Introduction 


A  useful  event  characterization  parameter  for  relatively  large  (e.g.,  mb  >  4)  teleseismic  events  is 
mb-Ms  (e.g.,  Marshall  and  Basham,  1972;  Blandford  et  al.,  1992),  provided  the  events  are  not  too 
deep.  To  use  mb-Ms  as  a  parameter  with  which  to  “screen”  events  of  natural  seismicity,  at  a  given 
confidence  level,  the  uncertainties  of  these  measurements  must  be  quantified.  Currently,  our 
treatment  of  the  uncertainty  in  mb-Ms  (e.g.,  Fisk  et  al.,  1995),  employed  by  event  screening  tools 
available  on  the  World  Wide  Web  site  at  the  Prototype  International  Data  Center,  Arlington,  VA, 
assumes  that  the  standard  deviations  for  individual  mb  and  Ms  measurements  are  0.3  magnitude 
units.  The  uncertainties  are  due  to  numerous  factors,  including  near-source,  path,  station, 
radiation  pattern  and  random  scattering  effects.  Many  of  these  factors  are  systematic  in  nature  and 
can  be  corrected  if  enough  information  is  available,  as  discussed  by  Rautian  and  Khalturin  (1994). 
Systematic  errors  due  to  radiation  pattern  effects  are  described  by  von  Seggem  (1970).1 

Our  current  approach  for  estimating  the  uncertainty  in  mb-Ms  further  depends  on  the  number  of 
observing  stations,  but  does  not  explicitly  account  for  correlations  between  station  measurements, 
which  are  due,  at  least  in  part,  to  radiation  pattern  effects.  In  this  report,  we  present  a 
generalization  to  our  current  approach  by  including  correlation  between  stations  due  to  station 
coverage.  Station  coverage,  to  be  defined  below,  is  essentially  how  much  azimuthal  angular 
coverage  the  stations  make  about  the  event  location.  This  effect  presumably  depends  on  the 
radiation  pattern  of  the  propagating  seismic  wave,  although  it  is  not  necessarily  the  case  that  such 
an  effect  can  be  seen  in  the  data.  The  methods  presented  below  to  treat  mb-Ms  uncertainty  due  to 
station  coverage  should  be  applied  in  computing  the  overall  mb-Ms  uncertainty  after  all  known 
corrections  for  systematic  errors  have  been  made.  A  similar  approach  that  treats  station 
correlation  for  mb  measurements  as  a  function  of  distance  on  the  focal  sphere  is  presented  in 
McLaughlin,  et  al.  (1988). 

In  the  remainder  of  Section  1,  we  describe  our  current  method  of  computing  a  confidence  interval 
for  mb-Ms.  In  Section  2,  we  describe  the  data  set  used  to  estimate  the  magnitude  uncertainties, 
based  on  all  events  in  the  Reviewed  Event  Bulletin  (REB)  for  which  Ms  was  measured  at  six  or 
more  stations.  In  Section  2,  we  also  describe  corrections  for  station  bias,  a  definition  of  station 
coverage,  and  estimates  of  the  standard  deviations  of  mb  and  Ms.  In  Section  3,  we  present  the 
multivariate  statistical  treatment  of  the  uncertainty  in  mb-Ms  due  to  station  coverage.  In  Section 
4,  we  provide  conclusions  and  recommendations  as  to  the  treatment  of  the  uncertainty  in  mb-Ms. 


1.  Such  corrections  are  being  pursued  by  the  Magnitude  Panel  (Murphy  et  al.),  established  by  the  Nuclear  Treaty 
Programs  Office.  We  plan  to  adopt  the  procedures  established  by  the  Magnitude  Panel  once  available. 
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1.1.  Current  Method  for  Computing  Magnitude  Uncertainty 


In  the  current  version  of  the  event  characterization  software,  installed  at  the  Center  for  Monitoring 
Research  (CMR),  an  event  retrieved  from  the  Reviewed  Event  Bulletin  (REB)  database  is  tested 
using  mb-Ms  in  the  following  way.  (A  test  for  focal  depth  is  also  applied  in  conjunction,  but  is  not 
discussed  here.)  For  an  event  for  which  there  are  both  mb  and  Ms  magnitude  measurements,  let 
Nb  be  the  number  of  stations  (nsta  in  the  schema)  with  mb  measurements.  The  network  average 
of  mb,  mb ,  is  given  by  the  sample  mean  of  the  Nb  data  values. 


Nh 


m 


b,j 


j=  1 


(1) 


where  the  sum  is  over  all  stations,  j ,  measuring  mb.  The  network  average  of  Ms,  Ms ,  is  computed 
in  a  slightly  more  complicated  manner  (Jepsen,  1996).  The  value  of  nsta  given  in  the  REB  is  the 
total  number  of  array  elements  measuring  Ms.  For  a  given  station,  j,  let  Uj  be  the  number  of 
elements  which  measured  Ms,  so  that 


Ns 

nsta  =  ^  Uj  ,  (2) 

j  =  1 

where  Ns  is  the  total  number  of  stations  measuring  Ms.  The  station  average,  Ms  is  given  by  the 
sample  mean  of  the  » .  values  of  Ms  at  station  j ,  and  the  network  average  of  Ms  is  then  given  by 
the  sample  mean  of  the  Ns  values  of  Ms  j. 


N. 


^  ^  X  ms,  j  • 


j  =  1 


(3) 


The  network  average  for  mb-Ms  is  then  taken  to  be  mb-Ms.  A  100(1  -a)%  confidence 
interval,  usually  taken  to  be  99%  (a  =  0.01),  about  the  expected  value  of  mb-Ms  is  then 
compared  with  1.2.  If  this  confidence  interval  is  entirely  below  1.2  (entirely  outside  the  nuclear 
explosion  population),  that  is,  if  the  largest  value  of  the  confidence  interval  is  less  than  1.2,  then 
the  event  is  considered  earthquakelike  and  screened.  The  confidence  interval  for  mb-Ms  is 
determined  as  follows. 


The  uncertainty  (standard  deviation)  for  both  individual  mb  and  Ms  measurements,  independent 
of  the  value  of  mb  or  Ms,  is  taken  to  be  0.3.  This  corresponds  to  a  variance  for  both  of  0.09,  that 

is. 


2 


(4) 


c\  =  c)  =  (0.3  )2  =  0.09. 

2 

The  variance  of  the  mean  of  Nb  independent  measurements  of  mb  is  given  by  Cb/Nb  with  a 

—  —  2 

similar  expression  for  Ms.  The  variance  of  mb-Ms,  denoted  by  c  ,  under  the  assumption  that 
mb  and  Ms  measurements  are  uncorrelated,  is 

c2  c2 

C2  =  variance^ -Ms)  =  jf  +  ■  (5) 

The  100 ( 1  -  a)%  confidence  interval  for  mb  -  Ms  is  determined  by  assuming  that  this  quantity 

2 

is  a  normally  distributed  random  variable  with  variance  a  with  a  true  (unknown)  mean  equal  to 
jx5.  Let  xa  be  defined  such  that  Prob(x>xa)  =  a  for  x  distributed  as  normal  with  zero 
mean  and  unit  variance.  Then 


Prob((ifc  -\is<mb-Ms  +  xac )  =  1  -  a  (6) 

and  the  one-sided  100(1  -a)%  confidence  interval  for  \lb  -  is  (-°° ,mb-Ms  +  xac]. 
Therefore,  an  event  with  measured  mb  -  Ms  is  screened  if 

mb-Ms  +  xac<  1.2  .  (7) 

(Using  a  two-sided  100(  1  -  a)%  confidence  interval,  [mb  -Ms- xa/2 c,  mb-Ms  +  xa/2c] ,  the 
screening  criterion  is  mb-Ms  +  xa/2 c <  1.2 .) 

Before  showing  how  this  magnitude  uncertainty  can  be  generalized  to  include  correlation  due  to 
station  coverage,  we  must  examine  the  REB  database  to  see  if  there  actually  is  such  correlation. 
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2.  Data  Analysis 

Before  describing  a  method  which  will  generalize  the  mb-Ms  magnitude  uncertainty  calculation 
it  is  necessary  to  examine  the  available  REB  data  to  determine  if,  indeed,  there  is  some 
dependence  in  the  variance  of  the  magnitudes  on  station  coverage.  An  appropriate  definition  of 
station  coverage  will  be  given  below. 

2.1.  Data 

For  the  purpose  of  studying  the  effect  of  station  coverage  on  magnitude  uncertainty  we  use  all 
events  in  the  REB  for  which  there  are  more  than  six  Ms  measurements  available.  There  were  400 
such  events,  dated  from  1 1  Jan  1995  (199501 1)  to  2  Jun  1996  (1996154).  Those  events  which  fell 
between  the  dates  17  Jan  1996  (1996017)  and  21  Feb  1996  (1996052)  were  not  used  due  to  an 
erroneous  formula  used  to  calculate  Ms  (Israelsson,  1996).  There  were  116  stations  in  this  list  of 
data.  Some  of  the  stations  consist  of  multiple  array  elements  each  of  which  may  measure  Ms  and 
these  individual  array  element  measurements  are  what  is  reported  in  the  stamag  table.  The  total 
number  of  array  elements  measuring  Ms  is  the  number,  nsta,  which  is  reported  in  the  netmag 
table,  and  not  the  total  number  of  stations  with  Ms  measurements.  The  network  average  Ms, 
reported  in  the  origin  and  netmag  tables,  average  over  stations  and  not  elements  and,  therefore,  is 
not  computed  by  summing  the  entries  in  the  stamag  table  and  dividing  by  nsta  (see  Eq  (3)  above 
for  the  correct  formula).  The  total  number  of  array  elements  in  the  list  is  141. 

For  each  event,  i,  the  sample  mean  for  mb  and  Ms  can  be  computed  as  in  Eq  (1)  and  (3).  The 
residual  for  a  given  mb  measurement  for  event  i  at  station  j  is  given  by 


eb,ij  =  mb,ij~mb,i > 


j  = 


b,  i  ’ 


(8) 


where  mh  ■■  is  the  mb  magnitude  for  event  i  at  station  j  and  ■  is  the  number  of  mb 
measurements  for  event  i .  The  residuals  for  Ms  measurements  are  given  in  a  similar  way.  The 
maximum  likelihood  estimate  (MLE)  of  the  variance  based  on  all  mb  measurements  is  given  by 


r 


,2  _ 


N 


b ,  i 


L  Vi  =  i 


-1  At. 


At, 


b,i 


I  I 


'b,ij 


J  J  «=  1  j= 1 


(9) 


with  a  similar  definition  for  sj .  Note  that  this  definition  is  simply  summing  the  residuals  for  all 
measurements  and  dividing  by  the  total  number  of  measurements,  which  is  sensible  only  if  it  is 
being  assumed  that  the  variance  of  mb  (and  Ms)  measurements  is  independent  of  magnitude 
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(independent  of  event).  When  Eq  (9)  (and  the  appropriate  equation  for  Ms)  are  applied  to  the  400 
selected  events  they  give 


sb  =  0.38,  ss  =  0.26.  (10) 

The  unbiased  sample  standard  deviations  for  mb  and  Ms,  which  are  obtained  by  subtracting  Ne  in 
the  denominator  of  Eq  (9),  are  sb  =  0.39  and  ss  =  0.28 . 

2.2.  Station  Corrections 

It  is  often  the  case  that  a  given  station  will  give  magnitude  readings  that  are  systematically  greater 
or  less  than  the  network  average  due  to  local  geophysical  properties  and  this  systematic  error 
accounts  for  some  of  the  variance  in  magnitude  measurements.  This  problem  is  currently  being 
addressed  by  Keith  McLaughlin  of  S-Cubed  (McLaughlin,  1996)  and  once  he  publishes  his  list  of 
station  corrections,  we  will  incorporate  them  into  the  Custom  Event  Characterization  Run.  For 
now,  we  will  estimate  station  bias  using  the  400  selected  events  from  the  REB.  We  define  the  bias 
at  station  j ,  bj ,  one  for  mb  and  one  for  Ms,  by 

bj=  (11) 

j  i 

where  the  sum  is  over  those  events,  total  number  equal  to  Nj ,  which  were  observed  at  station  j .  A 
corrected  magnitude  measure  at  station  j ,  one  for  mb  and  one  for  Ms,  is  then  given  by 

mj,  corrected  —  mj,  uncorrected  —  bj  *  (12) 

Using  corrected  values  for  mb  and  Ms  in  Eq  (9)  (and  the  appropriate  equation  for  Ms),  the 
variances  are  lowered  somewhat.  The  results  for  the  400  selected  events  are 

sb  =  0.33,  ss  =  0.23  (MLE);  sb  =  0.34,  ss  =  0.25  (unbiased  variance).  (13) 

2.3.  Definition  of  Station  Coverage 

Let  an  event  located  at  latitude  Xe  and  longitude  §e  be  observed,  by  mb  or  Ms,  at  N  stations  with 
latitudes  Xj  and  longitudes  (J )j,  j  =  1, ...,  N .  We  wish  to  define  a  number  q,  called  the  station 
coverage,  with  q  between  l/N  and  one,  such  that  q  near  one  corresponds  to  the  N  stations 
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being  in  directions  that  are  spread  uniformly  about  the  event  location  and  q  near  l/N 
corresponds  to  the  stations  all  in  the  same  direction.  It  is  easier  to  work  in  a  flat  x'y'  plane  than  on 
a  sphere  so  we  map  the  event  location  to  the  origin  of  the  x'y'  plane  and  the  N  station  locations  to 
points  on  the  unit  circle  about  the  origin,  with  station  one  on  the  x'  axis.  This  transformation, 
T:  (k,  <|))  ->  O',  /) ,  is  shown  in  Figure  1. 


/ 


Figure  1.  Mapping  of  event  and  station  locations  from  spherical  earth  to  x'y '  Cartesian  plane. 


To  carry  out  this  transformation  let  xe  be  a  unit  vector  in  the  direction  from  the  earth’s  center  to 
the  location  on  the  earth’s  surface  of  the  event.  In  a  right-handed  Cartesian  coordinate  system  with 
origin  at  the  earth’s  center,  z  axis  through  the  north  pole,  and  x  axis  through  zero  degrees 
longitude,  this  vector  is 

xe  =  cosA^cos^e^  +  cosA.gsin<])ee^  +  sinA,eez,  (14) 


where  ex ,  ey ,  and  ez  are  unit  vectors  along  the  x,  y  and  z  axes,  respectively.  Similarly, 

Xj  =  coskjCos§jex  +  coskjSin<\>jey  +  sin  Ayez ,  j  =  (15) 


is  a  unit  vector  in  the  direction  of  station  j .  A  plane  containing  the  origin,  the  event  location,  and 
station  j ,  has  a  unit  normal,  ny ,  perpendicular  to  both  xe  and  Xj ,  given  by 


x;xxe 


j  =  1, ...,  N . 


(16) 
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All  such  vectors  are  in  the  plane  normal  to  xe .  The  angle  made  by  two  great  circles  through  the 
event  location  and  two  stations  is,  by  definition,  the  angle  made  by  the  two  planes  containing  the 
origin,  the  event  location,  and  the  station  locations,  or,  equivalently,  the  angle  made  by  the 
normals  to  these  planes.  A  primed  coordinate  system  can  be  introduced  in  the  plane  normal  to  xe 
by  setting  the  unit  vector  along  the  x'  axis,  &x  ,  equal  to  fij ,  and  the  unit  vector  along  the  y'  axis, 
,  equal  to  a  vector  perpendicular  to  both  xe  and  &x  ,  that  is. 


e*  — 


n, 


=  xexex 


(17) 


In  this  coordinate  system,  the  points  with  coordinates 


j  =  1 . N, 


(18) 


are  points  on  the  unit  circle  with  point  1  on  the  x'  axis,  as  in  Figure  1. 

Consider  point  j  on  the  unit  circle  of  the  x'y'  plane,  which  makes  an  angle  a  •  with  the  x'  axis 
(0Cj  =  0 ).  The  sector  from  a;  -  p/2  to  (Xj  +  p/2 ,  as  shown  in  Figure  2,  has  total  angle  P . 


y' 


Figure  2.  Coverage  for  one  station,  equal  to  fi/27i  =  1/ N . 


The  fraction  of  the  total  angle,  A  =  2% ,  covered  by  the  angular  sector  of  size  P ,  is  p/A .  The 
interval  defined  by 
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(19) 


Ij  -  [a.j,bj]. 


where 


has  measure,  |l(7;)  =  (5/ A ,  which  is  less  than  or  equal  to  one  for  p  <  A . 

Now  consider  all  the  intervals,  Ij  =  [ap  bj],j=  1, . . N .  If  the  interval,  Ij ,  is  such  that  either 


dj  <  0  or  bj  >  1 ,  then  split  the  interval  into  two  intervals,  I  ■ 

=  7(^0  7(2),  where 

I j1)  =  [0,^],  /}2)  =  [1+^,1], 

if  a j<  0 , 

(21) 

/CD  =  [0  ,bj-l],  7(2)  =  [dj,  1] , 

if*;>l. 

(22) 

so  that  all  intervals  are  in  [0,  1  ] .  We  then  define  the  coverage,  q ,  of  N  stations  to  be 

<7  =  n(u/;l.  (23) 

S  =  i  J 

It  is  readily  seen  that  this  definition  yields  a  number  between  | 3/A  and  one.  If  (3  is  chosen  to  be 
A/N,  then  q  =  1  if  and  only  if  the  N  points  are  equally  spaced  on  the  unit  circle,  and 
q  =  1/N ,  the  smallest  possible  value,  if  and  only  if  all  stations  are  located  at  the  same  angle.  For 
example,  if  N  =  4  and  the  stations  are  located  at  0° ,  90° ,  180° ,  and  270° ,  then  q  =  1 ;  if  all 
the  stations  are  located  at  0° ,  then  q  =  0.25 ;  if  two  stations  are  at  0°  and  two  stations  are  at 
135° ,  then  q  =  0.50 ;  and  so  on. 

If  it  were  suspected  that  the  observed  radiation  pattern  for  the  magnitude  measurements  had  some 
symmetry,  then  it  might  be  useful  to  choose  A  to  be  something  other  than  360°  in  the  definition 
of  station  coverage.  For  example,  if  the  radiation  pattern  is  two-leafed,  i.e.,  symmetric  about  some 
axis,  then  the  magnitudes  at  points  180°  apart  would  be  the  same  (all  other  things  being  equal) 
and  so  the  coverage  might  more  usefully  be  defined  by  setting  A  =  180°,  which  considers 
intervals  separated  by  180°  to  be  the  same.  If  a  four-leafed  radiation  is  expected,  then  A  =  90° 
might  be  more  useful.  The  usefulness  of  choosing  A  to  be  something  other  than  360°  for  the 
actual  data  will  be  examined  below. 
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2.4.  Magnitude  Uncertainty  (uncorrected) 

To  study  the  effect  of  station  coverage  on  magnitude  uncertainty  we  examine  the  residuals  of 
subsets  of  the  data  (the  difference  between  the  mean  for  the  subset  and  the  mean  for  the  entire  set) 
as  a  function  of  station  coverage,  q.  For  a  given  event,  i,  let  Ni  be  the  number  of  station 
measurements  available  (for  either  mb  or  Ms)  and  let  mi  be  the  sample  mean  as  computed  by  Eq 
(1)  (for  mb)  or  Eq  (3)  (for  Ms).  Consider  all  subsets  of  the  Nt  stations  taken  n  at  a  time.  The  total 
number  of  such  subsets  is 


'N  a  =  AT,! 

n)  w!(iVf-n)!  * 


(24) 


For  example,  if  N(  =  5  and  n  =  2  then  there  are  10  subsets,  namely,  {1,2},  {1,3},  {1,4},  {1,5}, 
{2,3},  {2,4},  {2,5},  {3,4},  {3,5},  and  {4,5}.  For  a  given  subset,  {;'l5 ...,  jn} ,  denote  the  sample 
mean  relative  to  that  subset  by 


miUv-Jn)  =  -(miJi  +  ...+miJ),  l<jl<...<jn<Ni. 


(25) 


The  residual  for  each  subset  is  given  by 


•■■Jn)  =  ••• Jn)~mi ’  1  <  -<in^Ni- 


(26) 


Each  subset,  {jy, ...,  jn},  1  <  j1  <  ...  <jn  <  Nt,  consists  of  n  stations  surrounding  the  event 
location  with  a  corresponding  station  coverage,  q,  computed  using  Eq  (23).  In  studying  the 
relation  of  magnitude  uncertainty  and  station  coverage  it  would,  of  course,  be  better  to  use  the  true 
mean  rather  than  the  sample  mean  in  Eq  (26),  but  the  true  mean  is  not  known.  How  this  affects  the 
calculations  will  be  discussed  below. 

We  are  now  in  a  position  to  determine  if  there  is  any  dependence  of  uncertainty  on  station 
coverage.  Figure  3  is  a  scatterplot  of  the  square  root  of  the  square  (absolute  value)  of  the  residuals 
of  Ms  measurements  for  all  subsets  of  two  stations  for  every  event  in  the  database  satisfying  the 
conditions  that  the  number  of  Ms  stations  for  each  event  is  greater  than  six  and  the  station 
coverage  for  all  Ns  f  stations  used  to  compute  ms  ■  is  greater  than  0.6.  No  symmetry  is  assumed 
at  this  point  so  A  is  chosen  to  be  equal  to  360° . 
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Scatterplot  of  sdev(3/s)  vs.  q  (Ms),  (nsta(Ms)  =  2)\  (360°) 


q  (Ms) 


Figure  3.  Scatterplot  of  Ms  magnitude  uncertainty  versus  station  coverage  for  2  station  subsets,  A  -  360° . 


One  can  see  from  Figure  3  that  the  station  coverage  varies  from  q  =  1  /n  =  0.50  to  q  =  1 .  It  is 
difficult  to  tell  from  this  scatterplot  if  there  is  any  dependence  on  q.  To  aid  seeing  if  there  is  some 
dependence  on  q,  we  compute  the  average  variance  of  the  data  as  a  function  of  q.  To  do  this  we 
separate  the  data  into  K  bins,  qk  <  q  <  qk  +  j ,  k  =  1, K,  with  1/n  =  q^  <  ...  <  qK  + 1  =  1 . 
The  sample  variance  for  subsets  of  size  n  in  bin  k  is 


=  jr-  X  X  eUv-Jn)2’ 


'n,k  i=li<jl<...<jn<Ni 
?*<?<?*  + 1 


where  Ne  is  the  total  number  of  events,  the  second  sum  is  restricted  to  subsets  with  q  in  bin  k,  and 
Nn  k  is  the  total  number  of  terms  in  the  entire  summation.  Figure  4  contains  a  plot  of  s2>  k , 
k  =  1, ...,  5 ,  with  bins  [0.5, 0.6],  [0.6, 0.7],...,  [0.9,1. 0].  Each  snk  is  plotted  at  the  midpoint  of  its 
respective  bin.  The  error  bars  are  computed  by  assuming  the  k  are  distributed  as  chi-square 
with  Nn  k  degrees  of  freedom  and  choosing  a  95%  confidence  interval,  which  is  chosen  to  be 


[N n>  Icsn,  //Xo.025  ^n,  kSn,  k^0.915^  • 
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sdev(Afs)  vs.  q  (Ms),  (nsta(Ms)  =  2);  (360°) 


Figure  4.  Ms  magnitude  uncertainty  versus  station  coverage  for  2  station  subsets,  A  =  360°  (uncorrected). 


Figure  4  shows  a  dependence  on  q  which  has  a  maximum  at  q  -  0.50 ,  has  a  minimum  near 
q  =  0.75 ,  and  increases  a  q  approaches  one  but  not  to  the  maximum  value.  For  n  =  2  stations, 
q  =  0.50  corresponds  to  stations  at  the  same  angle,  q  -  0.75  corresponds  to  stations  separated 
by  90  degrees  and  q  =  1  corresponds  to  stations  separated  by  180  degrees.  It  is  seen,  then,  that 
the  greatest  uncertainty  occurs  when  two  stations  are  in  the  same  direction  from  the  event.  This 
implies  that  Ms  magnitude  measurements  from  stations  in  the  same  direction  are  positively 
correlated.  Minimum  uncertainty,  which  implies  a  negative  correlation,  occurs  when  stations  are 
separated  by  90  degrees.  Stations  separated  by  180  degrees  show  more  uncertainty,  more 
correlation,  than  those  separated  by  90  degrees  but  less  than  those  separated  by  0  degrees.  This 
suggests  that  there  may  be  a  radiation  pattern  with  a  180  degree  symmetry.  This  should  be 
observed  if  we  take  A  =  180°  in  our  definition  of  station  coverage,  q,  and,  as  Figure  5  shows, 
there  is  a  monotonic  decrease  of  magnitude  uncertainty  with  increasing  q. 

In  Figure  5,  q  =  0.50  corresponds  to  stations  separated  by  either  0  degrees  or  180  degrees,  while 
q  =  1  corresponds  to  separations  of  90  degrees.  If  there  were  a  prominent  four-leafed  radiation 
pattern  there  would  be  a  minimum  at  q  =  0.75 ,  corresponding  to  stations  separated  by  either  45 


ll 


degrees  or  135  degrees,  which  is  not  the  case.  A  plot  with  A  -  90° ,  which  is  shown  in  Figure  6, 
does  not  display  a  monotonic  decrease  in  uncertainty  with  q,  which  would  be  the  case  if  there 
were  prominent  four-leaf  radiation  patterns  present  in  the  Ms  magnitude  data. 


sdev(Ms)  vs.  q  (Ms),  (nsta (Ms)  =  2);  (180°) 


Figure  5.  Ms  magnitude  uncertainty  versus  station  coverage  for  2  station  subsets,  A  =  180°  (uncorrected). 


Before  discussing  the  relation  between  magnitude  uncertainty  and  station  coverage  for  more  than 
two  stations  and  for  mb  magnitude  data,  it  is  necessary  to  examine  the  effect  on  these  calculations 
of  using  finite-sized  data  samples,  and  making  any  necessary  corrections.  The  square  root  of  the 
variance  of  all  data  in  Figure  3  is  0.15.  Since  these  variances  were  computed  using  two  station 
combinations,  one  would  expect,  if  there  were  no  correlation  in  the  data,  for  the  variance  for  two 
stations  to  be  equal  to  the  variance  for  one  station  divided  by  two.  Or,  equivalently,  the  square  root 
of  two  times  0.15,  which  is  0.21,  should  be  equal  to  the  square  root  of  the  single  station  value,  ss , 
given  in  Eq  (13)  (we  are  using  data  corrected  for  station  bias,  here). 

This  value  is  0.23,  slightly  higher  than  0.21.  This  would  not  be  alarming  in  itself,  but  as  we  shall 
see,  as  the  number  of  elements  increase  in  the  stations  subsets,  the  difference  between  n  times  the 
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single  station  variance  and  the  computed  variance  increases  with  n,  which  is  due  to  the  use  of 
finite-sized  samples.  We  will  now  show  how  to  correct  the  variance  for  finite-sized  samples. 


sd ev(Ms)  vs.  q  (Ms),  (nsta (Ms)  =  2)\  (90°) 


Figure  6.  Ms  magnitude  uncertainty  versus  station  coverage  for  2  station  subsets,  A  =  90°  (uncorrected). 


2.5.  Correcting  for  Finite-Sized  Samples 

Consider  a  data  sample  of  size  N,  {xy,j  =  1, ...,  N} .  The  sample  mean,  as  usual,  is  given  by 

N 

(»> 

j=  1 

Define  the  sample  variance  over  all  subsets  of  size  n,  s% ,  by 


s 


2 

n 


»!(#-/»)! 

N\ 


+ 


1  <n<N. 


(30) 


For  n  =  1 ,  this  sum  reduces  to 
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(31) 


i  =  jj  Z<v-*>2  • 


;'=  i 

which  is  the  usual  maximum  likelihood  estimate  of  the  variance.  Eq  (30)  may  be  rewritten 

.2 


2  (W-/Q!  y 
ifl  N\  " 

ji*  •••*7„  =  iL 
N 

(N-n)\  1 
N\  2 


X ,  -X  X;  -X 

_i! - +  - 


n 


n 


X  [(^-  -  X)2  +  . . .  +  (xJn  -  X)2  +  2(X;.  -  x)(x;2  -*)  +  ...]. 
n  ji*  —  *jn  =  1 


Each  of  the  n  squared  terms  in  the  square  brackets  gives 

7  =  1 


ji*-*jn  =  1 
Eq  (32)  now  becomes 


N  N 

1  1  'V  /  — n  2  W  —  1  1  1 


(32) 


(33) 


7 1  ^  •••  *  In  ~  1 

and  each  of  the  n(n  -  1  )/2  remaining  cross  terms  gives 

jv  v 

^  2(Xjk-x)(xJm-x)  =  X  2(x--x)(xy.-x),  1  <k<m<n .  (34) 

^  1*7=1 


sn  =  n  N  2  +  I  (xi-WxrV  ■  (35> 

7=1  j  =  1 

To  evaluate  the  second  sum  in  Eq  (35)  note  that 


0  =  ( Nx-Nx )2  =  (xt  +  ...  +xN- Nx)2 


(36) 


giving 


=  (xj  -x+  ...  +xN-x) 

N  N 

=  X( xj~x)2+  X  (*,—  *)(*,—  *). 

y=l 


N 


N 


X  (*i  -  *) (*y  -  * )  =  -  X  -  • 

**7=1  7=1 


(37) 
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Finally,  we  have 


n  - 1 


Msifry-rf-fi- 


n  N-IJN 


j=  1 


n-  1 
N-l 


(38) 


Eq  (38)  shows  that  in  the  limit  of  N  large,  the  sample  variance  for  subsets  of  size  n  is  equal  to  the 
sample  variance  for  single  values  divided  by  n,  as  one  would  expected.  For  finite  N,  however, 
there  is  a  correction  factor,  less  than  one,  which  is  equal  to  zero  forn  =  N ,  as  it  should  be.  Using 
this  correction  factor  will  allow  us  to  get  a  better  estimate  of  the  dependence  of  magnitude 
uncertainty  on  station  coverage. 


2.6.  Magnitude  Uncertainty  (corrected) 

For  each  event,  i,  we  now  correct  the  variance  as  computed  in  Eq  (27),  by  using  Eq  (38)  to  adjust 
for  each  sample  of  size  Nt  used  to  compute  mi .  The  corrected  sample  variance  for  subsets  of  size 
n  in  bin  k  is  then 


k\  corrected 


cik<<l<clk+  1 


(39) 


Figure  7  is  a  plot  of  s2>  k;  corrected  similar  to  Figure  4.  Notice  that  the  plot  is  basically  shifted  up  by 
a  small  amount.  In  fact,  the  average  square  root  of  the  variance  of  this  curve  times  n  (=  2)  is  0.23, 
now  equal  to  the  single  station  value,  ss ,  given  in  Eq  (13). 

Figures  8,  9,  and  10  plot  sn  ^corrected  f°r  subsets  with  n  =  3, 4,  5 ,  respectively,  and  A  =  360° , 
that  is,  no  symmetry  is  assumed.  For  given  n,  only  those  events  which  have  Ns  greater  than  n+  4 
are  used,  which  reduces  the  number  of  events  used  as  n  increases.  In  each  case,  it  can  be  seen  that 
the  general  trend  is  for  a  decrease  in  Ms  magnitude  uncertainty  as  a  function  of  increasing  q.  This 
is  the  expected  result;  the  more  the  stations  are  evenly  distributed  around  the  unit  circle,  greater  q, 
the  less  the  expected  variance  in  magnitude.  The  average  variance  for  each  case  is  approximately 
equal  to  the  single  station  value  divided  by  n,  as  expected  from  Eq  (38).  Note  that  the  large  error 
bars  for  the  endpoints  on  some  of  the  plots  are  due  to  the  small  number  values  which  fall  into  the 
particular  bin  used  to  calculate  the  average.  The  seeming  great  decrease  in  uncertainty  as  a 
function  of  q  for  n  =  5  should  not  be  taken  too  seriously  due  to  the  large  error  bars  on  the 
endpoints. 
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sd ev(Ms)  vs.  q  (Ms),  (nsta (Ms)  =  2);  (360°) 


Figure  7.  Ms  magnitude  uncertainty  versus  station  coverage  for  2  station  subsets,  A  =  360°  (corrected). 


q  (Ms) 

Figure  8.  Ms  magnitude  uncertainty  versus  station  coverage  for  3  station  subsets,  A  =  360°  (corrected). 
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s dev(Ms)  vs.  q  (Ms), 


Figure  9. 


Figure  10. 


(nsta(Ms)  =  4)\  (360°) 


Ms  magnitude  uncertainty  versus  station  coverage  for  4  station  subsets,  A  =  360°  (< 


sd ev(Ms)  vs.  q  (Ms),  (nsta (Ms)  =  5);  (360°) 


q  (Ms) 

Ms  magnitude  uncertainty  versus  station  coverage  for  5  station  subsets,  A  =  360°  (< 
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sdev(.Ms)*sqrt(n)  vs.  q  (Ms),  (nsta (Ms)  =  1 ,2,3,4, 5);  (360°) 


Figure  11.  Ms  magnitude  uncertainty  times  the  square  root  of  number  of  stations  versus  station  coverage  for 
n  =  1, 2, 3.4,  5  ,A  =  360°  (corrected). 

Figure  11  is  a  plot  of  Jn  x  snk.c orrected  as  a  function  q  for  subsets  with  n  =  1, 2,  3, 4,  5  and 
A  =  360° ,  where  the  error  bars  have  been  suppressed  for  clarity.  This  plot  combines  each  of  the 
plots  in  Figures  7-10  times  the  square  root  of  n.  As  expected,  each  of  the  curves  has  about  the 
same  average  variance,  approximately  equal  to  ss  =  0.23 ,  which  is  displayed  as  a  horizontal  line 
in  the  plot  and  labeled  n  =  1 .  Recalling  that  the  endpoints  were  computed  with  very  few  data 
points  and,  therefore,  have  large  error  bars  associated  with  them,  the  curves  have  approximately 
the  same  dependence  on  q,  which  implies  that  we  have  chosen  a  consistent  definition  of  station 
coverage.  In  review,  recall  that  these  curves  were  computed  by  using  Ms  data  from  400  events  in 
the  REB,  which  have  been  corrected  for  station  bias.  The  bias  was  estimated  from  these  events. 
The  sample  variance  for  subsets  of  a  given  size,  n,  were  corrected  for  finite  sample  size,  as  in  Eq 
(39),  and  each  curve  was  then  adjusted  by  multiplying  by  the  square  root  of  n.  What  remains 
should  then  be  a  good  estimate  of  the  true  dependence  of  the  Ms  magnitude  uncertainty  on  station 
coverage.  If  the  definition  of  station  coverage  is  chosen  well,  the  station  coverage  dependence 
should  be  independent  of  n,  as  is  approximately  the  case  in  Figure  1 1.  In  Figure  12  we  combine 
Figures  7-10  for  comparison,  which  have  been  off-set  for  clarity. 
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sdev(jl/s)  vs.  q  (Ms),  (nstafAfsj  =  1 ,2, 3,4,5);  (360°) 


Figure  12.  Ms  magnitude  uncertainty  versus  station  coverage  for  n  =  1,2, 3.4,  5 ,  A  =  360°  (corrected). 


Figure  13  plots  sn  corrected  as  a  function  of  q  for  subsets  with  n  =  1,2,  3, 4,  5  and  A  =  180° , 
which  assumes  that  the  data  is  symmetric  with  respect  to  180  degree  rotations.  It  can  be  seen  that 
for  each  n  greater  than  one  there  is  a  general  decrease  in  Ms  magnitude  uncertainty  with 
increasing  q,  but  that,  compared  to  Figure  12,  the  slopes  are  not  as  steep.  Although  this  is  mild 
evidence  that  there  is  a  180  degree  symmetry  in  the  Ms  data,  it  seems  preferable  not  to  make  this 
assumption  when  attempting  reduce  Ms  magnitude  uncertainty  when  stations  are  well  spread 
about  the  event  location. 

mb  Magnitude  Uncertainty 

We  now  consider  the  relation  between  mb  magnitude  uncertainty  and  station  coverage  in  an 
analogous  manner  as  the  Ms  analysis  described  in  detail  above.  Figures  14  and  15  are  plots  of 
sn  £•  corrected  for  mb  measurements  as  a  function  of  q  for  subsets  with  n  =  1,2,  3, 4,  5  for 
A  =  360°  and  A  =  180°,  respectively.  Unlike  the  Ms  case  shown  in  Figure  12,  the  mb 
magnitude  uncertainty  for  A  =  360°  does  not  show  consistent  decreasing  trends  as  q  increases, 
although  there  seems  to  be  slight  indications  of  decreasing  trends  when  a  180  degree  symmetry  is 
assumed  as  shown  in  Figure  15. 
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sdev(Ms)  vs.  q  (Ms),  (nsta (Ms)  =  1 ,2,3,4, 5);  (180°) 


Figure  13.  Ms  magnitude  uncertainty  versus  station  coverage  for  n  =  1, 2, 3.4,  5 ,  A  =  180°  (corrected). 


It  should  be  noted  that  if  the  dependence  on  station  coverage  is  due  to  the  radiation  pattern  of  the 
seismic  wave,  and  this  is  not  necessarily  the  case,  then  for  mb  magnitudes,  which  are  produced  by 
body  waves,  such  a  pattern  would  most  readily  be  observed  by  mapping  station  locations  to  the 
focal  sphere  (see,  for  example,  Aki  and  Richards,  1980).  It  is  not  clear  that  large  q  for  mb 
measurements  on  the  earth’s  surface  correspond  to  large  coverage  of  the  focal  sphere.  To 
determine  this,  one  would  have  to  use  raytracing  techniques,  which  are  beyond  the  scope  of  this 
study.  The  mild  trends  seen  in  the  A  =  180°  case  are  not  completely  surprising,  however, 
because  any  180  degree  symmetry  in  the  radiation  pattern  on  the  focal  sphere  would 
approximately  propagate  to  the  surface  due  to  the  approximate  spherical  symmetry  of  the  earth’s 
interior.  In  any  case,  the  dependence  on  q  is  small  enough  that  for  practical  purposes  we  will  say 
that  there  is  no  dependence  of  mb  magnitude  uncertainty  on  surface  station  coverage. 
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sdev(mb) 


0.4 


sdev(mb)  vs.  q  (mb),  (nsta(mb)  =  1 ,2,3,4, 5);  (360°) 


Figure  14.  mb  magnitude  uncertainty  versus  station  coverage  for  n  =  1,  2, 3.4, 5 ,  A  =  360°  (corrected). 


sdev(m6)  vs.  q  (mb),  (nstafmbj  =  1 ,2,3,4>5)\  (180°) 


Figure  15.  mb  magnitude  uncertainty  versus  station  coverage  for  n  =  1, 2, 3.4, 5 ,  A  =  180°  (corrected). 


3.  Multivariate  Approach 

In  the  previous  section  it  was  shown  how  the  variance  in  Ms  and  mb  measurements  depends  on 
station  coverage.  Generally  speaking,  as  the  stations  observing  Ms  are  more  uniformly  spread 
about  the  unit  circle  the  smaller  the  variance.  The  dependence  of  the  variance  of  mb 
measurements  on  station  coverage  is  substantially  less  than  for  Ms  measurements.  Inspection  of 
Figure  11,  which  plots  (corrected)  Ms  magnitude  uncertainty  (standard  deviation)  times  the 
square  root  of  n,  the  number  of  stations,  as  function  of  q,  the  station  coverage,  for  n  equal  to  1,2, 
3, 4,  and  5,  shows  that  the  dependence  of  the  Ms  uncertainty  on  q,  once  scaled  by  the  square  root 
of  n,  is  essentially  independent  of  n.  This  suggests  that  the  portion  of  the  variance  of  the  mean  of 
Ms  for  a  given  set  of  n  measurements,  which  depends  on  station  coverage,  can  be  determined  by 
the  relative  correlation  between  each  pair  of  stations,  as  a  function  only  of  the  angle  made  by  the 
two  stations  with  respect  to  the  event  location.  Whether  this  assumption  is  warranted  or  not  for  the 
actual  REB  data  will  be  addressed  below.  A  multivariate  approach  is  most  suited  in  this  case. 

3.1.  Uncertainty  in  mb-Ms  Using  Multivariate  Statistical  Analysis 

We  wish  to  generalize  the  expression  for  the  variance  of  mb  -  Ms  given  in  Eq  (5)  by  considering 
the  correlations  between  stations  for  Ms  measurements  as  a  function  of  station  angle.  We  will 
assume  that  there  is  no  correlation  in  the  individual  mb  measurements,  as  discussed  above,  and 
that  there  is  no  correlation  between  Ms  and  mb  measurements.  In  fact,  for  the  400  events  we  are 
considering  in  the  REB,  the  sample  correlation  coefficient  between  Ms  and  mb  measurements, 
using  those  station  measurements  that  have  both  Ms  and  mb  readings,  is  approximately  0.06,  so 
this  assumption  is  well  justified.  Under  these  assumptions,  the  variance  of  mb-Ms,  denoted  by 
a2 ,  is  given  by 

2 

2  _  (7  l  .  _ 

a  =  variance  (mfc  -  Ms)  =  —  +  variance(MJ) .  (40) 

Nb 

To  calculate  the  variance  of  Ms ,  suppose  that  for  a  given  event  there  are  n  Ms  measurements 
assumed  to  be  given  by  the  n  components  of  a  random  vector  x  =  (xj , . . . ,  xn)' ,  such  that 

Xj  =  m  +  ej,  j  =  1, n ,  (41) 

where  m  is  the  mean  of  each  measurement,  assumed  to  be  the  same  for  all  j  =  1, ...,  n ,  and 
E(ejek)  =  'Zjf. .  In  vector  notation  we  have 
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x  =  ma  +  e, 


E(ee’)  =  Z, 


(42) 


where  a  is  a  vector  with  all  components  equal  to  one,  a  =  ( 1, . . 1 )' .  Any  linear  combination  of 
the  Xj  of  the  form 


z  -  c'x,  with  c'a  =  1, 

(43) 

has  mean  equal  to  m,  E(z)  =  m,  and  variance  equal  to  c'Zc .  In  particular,  if  c  = 
is  equal  to  the  usual  sample  mean,  Ms ,  with  variance 

(l/n)a,  thenz 

variance  (A/)  =  a'Za , 

Nj 

(44) 

where  we  have  now  written  Ns  for  n. 

We  also  note  for  reference  that  if  c  =  (a'Z  1  a)  1 Z  1  a ,  then  z,  under  the  assumption  of  normality, 

gives  the  maximum  likelihood  estimate  (MLE)  of  m,  which  is  also  the  linear  combination  having 
~  .  _1  -1 
mean  m  with  minimum  variance,  equal  to  (a'Z  a)  . 

We  now  assume  that  for  a  given  event,  the  covariance  matrix,  Z ,  has  the  form 

^jj  =  °s  ’  J  =  ~'’NS> 

(45) 

^jk  =  Pjk^h  j*k  =  1,  ...,NS, 

(46) 

where 

pjk  =  f(cos(ak-a.j)),  (47) 

and  a  •  is  the  azimuthal  angle  of  station  j  relative  to  the  event  location.  That  is,  the  variance  of 
each  Ms  measurement  at  each  station  is  assumed  to  be  the  same,  whose  value  is  given  in  Eq  (10) 
when  corrections  for  station  bias  are  not  used  or  Eq  (13)  if  corrected  for  bias,  and  the  correlation 
coefficient,  p  -k ,  depends  only  on  the  angular  separation  of  station  j  and  k  with  respect  to  the  event 
location.  The  function /will  be  estimated  from  the  REB  data  in  the  next  section. 
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3.2.  The  Correlation  Coefficient 


To  estimate  the  function /in  Eq  (47)  we  can  use  the  curve  plotted  in  Figure  7  for  Ms  magnitude 
standard  deviation  as  a  function  of  q  for  two  station  subsets.  From  the  definition  of  q  for  two 
stations  it  is  readily  seen  that  a  =  2it (q  -  1/2) ,  where  a  is  the  angle  made  by  the  two  stations 
with  respect  to  the  event  location.  To  convert  the  variances  in  Figure  7  to  correlation  coefficients 
note  that  for  the  random  variable  ( x  +  y)/ 2 , 

+  >0/2  =  ax/2  +  tf/2  +  2c0V(f  >  =  \*l  +  \°y  +  f C0V(*>  >)  *  (48) 

Assuming  that  o%  =  <Jy  =  of,  we  have 

=  COV(x,  y)  _  2  a(x  +  y)/2  _  1  (49) 

Clj 

Using  this  transformation  we  plot  in  Figure  16  the  correlation  coefficient  as  a  function  of  cosine 
of  the  angle  between  two  stations  for  all  two  station  subsets  in  the  REB  database  with  total 
number  of  Ms  measurements  greater  than  six.  The  data  has  been  split  into  five  bins. 

To  select  the  function/,  we  fit  the  five  points  in  Figure  16  with  a  function  of  the  form 

2 

/( cosa)  =  b0  +  bjcosa  +  b2cos  a.  (50) 


Using  least  squares  the  coefficients  are 

b0  =  -0.17,  b1  =  0.13,  b2  =  0.35 .  (51) 

The  function  /(cosa)  is  also  plotted  (the  smooth  curve)  in  Figure  16.  Notice  that  for  small 
station  angles,  cosine  near  one,  the  correlation  is  positive,  suggesting  that  stations  which  have 
angles  from  an  event  which  are  close  have  Ms  measurements  that  tend  to  be  above  average  or 
below  average.  For  station  angles  close  to  90  degrees,  cosine  near  zero,  the  correlation  coefficient 
is  negative,  consistent  with  radiation  patterns  that  are  larger  than  average  in  a  given  direction  and 
smaller  than  average  in  a  direction  90  degrees  from  the  first  direction.  The  correlation  coefficient 
for  stations  diametrically  opposed,  cosine  near  minus  one,  is  nearly  zero.  If  the  effect  of  radiation 
patterns  with  180  degree  rotational  symmetry  were  strong,  we  would  expect  this  value  to  be 
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positive,  the  same  as  zero  degrees.  Since  this  is  not  the  case,  we  will  make  no  assumptions  about 
the  symmetry  properties  of  the  correlation  coefficient  function. 


Correlation  vs.  Cosine  of  Station  Angle 


Figure  16.  Correlation  coefficient  as  a  function  of  cosine  of  the  station  angle  for  2  station  subsets  and  smooth 
fit 


3.3.  The  Consistency  of  the  Correlation  Coefficient  Function 

It  will  be  recalled  that  we  have  assumed  that  the  correlation  of  two  Ms  magnitude  measurements 
for  a  single  event  is  independent  of  event  location  and  depends  only  on  the  angle  the  two  stations 
make  with  the  event  location.  If  this  were  not  true,  then  it  would  be  possible  for  the  function  given 
in  Eq  (50),  which  was  estimated  using  actual  data,  to  give  inconsistent  results,  that  is,  it  would  be 
possible  for  the  covariance  matrix,  E ,  defined  by  Eq  (45)  and  (46)  to  not  be  positive  definite, 
which  it  must  be  if  it  were  a  true  covariance  matrix.  If  we  can  prove  that  E  is  always  positive 
definite  for  all  choices  of  station  angles,  then  our  correlation  function  is  consistent,  although,  of 
course,  it  does  not  prove  that  the  assumption  that  the  correlation  depends  only  on  station  angle  is 
true.  We  will  prove  below  the  somewhat  weaker  result,  although  it  may  be  equivalent,  that  the 
variance  of  Ms ,  as  given  in  Eq  (44),  is  always  positive. 
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Using  a  =  (1, 1)'  Eq  (44)  gives 


variance^)  =  -^a'Ea  =  ^[l+^f  X  Pj*l  • 

sL  s^i<j<k<Ns  JA 

Defining 

F(av ...,  an)  =  X  P jk=  X  /(cos(cc*-ay.)), 

1 <j<k<n  1 <j<k<n 

we  must  show  that  F(av  ...,  an)  >  -n/2  for  all  a^,  j  =  1, n  and  all  n. 
We  first  show  that  if  fix)  is  convex  for  -1  <x<  1 ,  that  is, 

fiXx  +  ( 1  -  X)y)  <  ff(x)  +  ( 1  -  X )f(y) , 


(52) 


(53) 


(54) 


whenever  -1  <x<y<  1  and  0  <  X  <  1 ,  then  F  achieves  its  minimum  when  the  a ,  •  are  evenly 
spaced.  For  functions  with  continuous  second  derivatives,  convexity  is  equivalent  to  having 
f"(x)  >0  for  -1  <  x  <  1 .  Since,  as  is  easily  shown,  F(aj  +  a, ...,  an  +  a)  =  F(av  a„)  for 
all  a,  it  is  sufficient  to  show  that  F(0, 2%/n, ...,  (n-  l)27t/n)  <F(al5 ...,  a„)  for  all  ay-.  To 
see  this,  notice  that  if  some  a  j ,  say  ,  was  not  equal  to  0  then  half  of  the  terms  in  Eq  (53)  which 
contain  will  have  arguments  of/ which  increase,  by  Ax,  say,  and  the  other  half  will  have 
arguments  that  decrease  by  Ax ,  because  the  points  are  evenly  spaced  on  the  unit  circle  (if  n  is 
even  the  point  that  is  diametrically  opposed  can  be  written  /  =  f/2  +  f/2  and  the  argument  is 
the  same).  By  the  convex  property  of/ with  X  =  1/2 , 

f{x  -  Ax)  +  f{x  +  Ax)  >  2/ (x) ,  (55) 


which  shows 

F(0, ...,  ij+\)2%/n, ...,  in-  \)2%/n)  < F(0, ...,  (;  +  1  )2n/n  +  a, ...,  in-\)2n/n)  (56) 

for  all  a .  The  extension  to  two  variables  at  a  time  and  finally  to  all  n  variables  is  similar  but 
tedious.  For  example,  for  two  variables,  inequalities  such  as 
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f(x  +  Axt )  +  f(x  -  Ax{  +  Ax2)  +  f(x  -  Ax2)  ^  2  f(x  +  Ax2/2)  +  f(x  -  Ax2 )  (57) 

>  2  f(x  +  Ax2/2)  +  2  f(x  -  Ax2/2)  -  f(x)  >  3  f(x), 

must  be  employed. 


Having  shown  that  F  achieves  its  minimum  when  the  a ,  •  are  evenly  spaced,  we  must  have 

W  “•  1 

F(0,  2%/ n,  12tt/ n)  =  X  (n -  j)/jcos|y'^J)  >-| .  (58) 

j-  1 


That  is,  the  convex  function,/,  must  satisfy 

n 
2 


2-i 


X  (cos(/vI) + *  for  n  even> 

7=1 


and 


«-i 

2 

X  /(cos  (7^  J) >  ~\ '  for  n  odd- 

7=1 


(59) 


(60) 


Note,  these  equations  show  that /can  not  be  entirely  negative,  i.e.,  Ms  measurements  can  not  be 
negatively  correlated  for  all  angles. 


It  now  remains  to  show  for/ given  in  Eq  (50)  and  (51),  which  is  clearly  convex,  Eq  (59)  and  (60) 
are  satisfied  for  all  n.  Figure  17  plots  the  left  hand  sides  of  Eq  (59)  (for  n  even)  and  Eq  (60)  (for  n 
odd)  for  n  -  1, ...,  30 .  It  is  apparent  that  sums  are  well  above  -1/2  for  all  n.  To  show  this  for 
large  n,  denote  the  left  hand  side  of  Eq  (60)  by  sn .  Then 

b0  +  b1cos(j^j  +  b2cos  (j^y  (61) 

For  large  n, 

n/2  2  n 

X  cos  ^ j~~j  ~  ~  J  cos  2xdx  =  ^ ,  (62) 

7=1  0 
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3.4.  Results 


We  now  apply  our  formula  for  the  variance  of  mb-Ms,  given  in  Eq  (40)  and  using  Eq  (44)-(47), 
(50)  and  (51),  to  all  400  events  in  the  REB  database  with  the  number  of  Ms  measurements  greater 
than  six.  For  each  event  we  calculate  the  variance  of  Ms  as  in  Eq  (44)  and  the  variance  of 
mb-Ms  as  in  Eq  (40).  The  square  root  of  the  average  variance  of  Ms  times  Ns  over  the  400 
events  was  0.34.  If  there  had  been  no  correlation  between  stations,  this  value  would  have  been 
0.28,  the  unbiased  sample  standard  deviation.  We  can  see  that  the  value  with  correlation  is  larger, 
indicating  that  on  average,  the  events  in  the  data  base  have  positive  correlation.  The  square  root  of 
the  average  variance  of  mb  times  Nb  plus  the  average  variance  of  Ms  times  Ns  over  the  400 
events  was  0.52,  compared  to  0.48  if  there  had  been  no  correlation.  Thus,  on  average  the  effect  of 
correlation  is  small,  but  can  vary  for  different  events,  depending  on  the  azimuthal  coverage  of  the 
stations  which  recorded  them.  If  Ns  =  4 ,  for  example,  the  minimum  value  of  jNsvai(Ms)  is 
0.84  •  Gs  and  the  maximum  value  is  1.39  ■  Gs ,  a  range  in  values  of  55%. 
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4.  Conclusions  and  Recommendations 


Until  the  magnitude  panel  has  provided  its  conclusions  on  magnitude  estimates  and  their 

uncertainties,  we  recommend  using  the  following  formulas  to  estimate  the  uncertainty  in  mb-Ms. 
_  2 

The  variance  of  mb  -  Ms ,  denoted  by  a  ,  is  given  by 


a =  +  variance(Mj) , 


(65) 


where 


variance 


iViL  1  s\<j<k<Ns  /J 


/(cosa)  =  -0.17  +  0.13cosa  +  0.35cos  a. 


(66) 

(67) 


Here,  Nb  is  the  number  of  mb  measurements,  Ns  is  the  number  of  Ms  (station)  measurements, 
and  a .,  j  =  1, Ns,  is  the  azimuthal  angle  with  respect  to  the  event  location  of  Ms 
measurement  j.  If  station  measurements  are  not  corrected  for  bias,  we  recommend 

Gb  =  0.39,  cs  =  0.28.  (68) 

If  each  station  measurement  is  corrected  by  the  average  bias  for  current  data  in  the  REB,  we 
estimate  the  corrections  will  give 

Gb  =  0.34,  gs  =  0.25.  (69) 

We  anticipate  that  these  values  will  be  improved  by  the  magnitude  panel  with  all  appropriate 
systematic  errors  removed. 

Finally,  an  event  with  measured  mb  -  Ms  is  screened  at  the  100(  1  -  a)%  confidence  level  if 

mb-Ms  +  xaG<  1.2,  (70) 

where  xa  is  defined  such  that  Prob(r  >xa)  =  a  for  x  distributed  as  normal  with  zero  mean  and 
unit  variance. 
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